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The value of the V u d matrix element of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix can be derived from nuclear superallowed beta decays, neutron decay, and 
pion beta decay. We survey current world data for all three. Today, the most 
precise value of V u d comes from the nuclear decays; however, the precision is 
limited not by experimental error but by the estimated uncertainty in theoretical 
corrections. The neutron data arc approximately a factor of four poorer in precision 
but this could change dramatically in the near future as planned experiments 
come to fruition. The nuclear result (and the most recent of the neutron decay 
results) differ at the 98% confidence level from the unitarity condition for the CKM 
matrix. We examine the reliability of the small calculated corrections that have 
been applied to the data, and assess the likelihood of even higher quality nuclear 
data becoming available to confirm or deny the discrepancy. Some of the required 
experiments depend upon the availability of intense radioactive beams. Others are 
possible today. 



1 Introduction 

The Cabibbo-Kobayashi-Maskawa matrix 0i relates the quark eigenstates of 
the weak interaction with the quark mass eigenstates (unprimed) 

(Vud Vus V ub \ fd\ 

= Vcd V cs V cb \ [ s\ (1) 

V Vtd Vts VtbJ \bj 

and, as such, the matrix is unitary. Thus there are many relationships among 
the nine elements of the matrix that can be tested by experiment. The leading 
element, V ud , only depends on quarks in the first generation and so is the 
element that can be determined most precisely. Here we will discuss the current 
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status of V uc i and the unitarity test as it relates to the elements in the first 
row: 

v& + v£ + v& = i. (2) 

In examining the unitarity test, we will adopt the Particle Data Group 
(PDG98)t3 recommendations for V us and V u b- We note in particular that, in 
its 1998 update, PDG98 is recommending for V us only the value determined 
from K e 3 decay, | V us |= 0.2196 ± 0.0023, arguing that the value obtained 
from hypcron decays suffers from theoretical uncertainties due to first-order 
SU(3) symmetry-breaking-jeffects in the axial-vector couplings. In his talk to 
this conference, MarcianoQ notes that new experimental studies of K e ^ decay 
are underway (existing data are 20 years old) and SU(3) symmetry-breaking 
effects, which in this case are of second order, are being reexamined. As to V u b, 
its value is small, | V u b |= 0.0032 ± 0.0008, and consequently it has a negligible 
impact on the unitarity test, Eq. (||). 

2 The value of V ud 

The value of V u a can be determined from three distinct sources: nuclear super- 
allowed Fermi beta decays, the decay of the free neutron, and pion beta decay. 
We discuss each in turn. 

2.1 Nuclear superallowed Fermi beta decays 

Nuclei have the singular advantage that transitions with specific characteristics 
can be selected and then isolated for study. One example is the superallowed 
+ — > + beta transitions, which depend uniquely on the vector part of the 
weak interaction. Furthermore, in the allowed approximation, the nuclear 
matrix element for these transitions is given by the expectation value of the 
isospin ladder operator, which is independent of any details of nuclear struc- 
ture and is given simply as an SU(2) Clebsch-Gordan coefficient. Thus, the 
experimentally determined /i-values are expected to be very nearly the same 
for all + — > + transitions between states of a particular isospin, regardless 
of the nuclei involved. Naturally, there are corrections to this simple picture 
coming from electromagnetic effects, but these corrections are small - of or- 
der 1% - and calculable. Thus, if we write 5r as the nucleus-dependent part 
of the radiative correction, as the nucleus-independent part of the radia- 
tive correction, and 5c as the isospin symmetry-breaking correction, then the 
experimental /lvalue can be expressed as follows: 
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Table 1: Experimental results (Qec> ti/2 an< ^ branching ratio, R) and calculated correction, 
PeCj f° r + — > 0+ transitions. The other calculated corrections, <5fj and 8q, are given in 
Tables H and W respectively. 





Qec 


h/2 


R 


Pec 


ft 






(keV) 


(ms) 


(%) 


(%) 


w 


M 


10 c 


1907.77(9) 


19290(12) 


1.4645(19) 


0.296 


3038.7(45) 


3072.9(48) 


14 Q 


2830.51(22) 


70603(18) 


99.336(10) 


0.087 


3038.1(18) 


3069.7(26) 


26m js^ 


4232.42(35) 


6344.9(19) 


> 99.97 


0.083 


3035.8(17) 


3070.0(21) 


34 C1 


5491.71(22) 


1525.76(88) 


> 99.988 


0.078 


3048.4(19) 


3070.1(24) 


38m K 


6044.34(12) 


923.95(64) 


> 99.998 


0.082 


3049.5(21) 


3071.1(27) 


42 Sc 


6425.58(28) 


680.72(26) 


99.9941(14) 


0.095 


3045.1(14) 


3077.3(23) 


46 V 


7050.63(69) 


422.51(11) 


99.9848(13) 


0.096 


3044.6(18) 


3074.4(27) 


so Mn 


7632.39(28) 


283.25(14) 


99.942(3) 


0.100 


3043.7(16) 


3073.8(27) 


54 Co 


8242.56(28) 


193.270(63) 


99.9955(6) 


0.104 


3045.8(11) 


3072.2(27) 










Average, Tt 


3072.3(9) 












x 2 /» 


1.10 



ft(l + 8 R )(1 - S c ) =Tt= 2G 2 (1 + AV) = constant, (3) 

where G v is the weak vector coupling constant and K = 2ir 3 In 2fi(tic) 6 / (m e c 2 ) 5 , 
which has the value K/(hc) 6 = (8120.271 ± 0.012) x 10~ 10 GeV- 4 s. Thus, to 
extract V u d from experimental data, the procedure is to determine the Tt- 
values for a variety of different nuclei having the same isospin, and then to 
test if they are self-consistent. If they are, their average is used to determine 
a value for G v and, from it, V u d- 

To date, superallowed + — * + transitions have been measured to ±0.1% 
precision or better in the decays of nine nuclei ranging from 10 C to Co. 
World data on Q-values, lifetimes and branching ratios - the results of over 
100 independent measurements - were thoroughly surveyed^ in 1989 and then 
updated again El for the last WEIN conference. Since that time, there has 
been a new n measurement for the 10 C branching ratio, which came from an 
experimentQ using Gammasphere, and one forlhe Q-value of the 38rra K decay, 
a result reported to this conference by Barker a . We have incorporated both 
these new measurements into our data base, and list the resulting weighted 
averages in the first three columns of Table [l| Using the calculated electron- 
capture probabilities □ given in the next column, we obtain the "uncorrected" 
ft- values listed in column 5 with partial half-lives determined from the formula 
t = t 1/2 (l + P EC )/R. 
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Figure 1: Tt- values for the nine precision data, and the best least-squares one-parameter fit. 



We save a detailed discussion of the radiative and Coulomb corrections for 
Sec. H In the present context, though, it should be noted that the values for 
5r and 5c used to derive Tt were taken from the last column of Tables |3| and 
|], respectively; each value results from more than one independent calculation. 
In the case of 5r, the calculations are in complete accord with one another; 
for 5c, we have used an average of two independent calculations with assigned 
uncertainties that reflect the (small) scatter between them. Thus, in a real 
sense, both experimentally and theoretically, the ^-values given in Table ^ 
and plotted in Fig. [I] represent the totality of current world knowledge. The 
uncertainties reflect the experimental uncertainties and an estimate of the rela- 
tive theoretical uncertainties in 5c There is no statistically significant evidence 
of inconsistencies in the data (x 2 / v — 1-1), thus verifying the expectation of 
CVC at the level of 3 x 10 -4 , the fractional uncertainty quoted on the average 
Tt- value. 

In using this average Tt- value to determine G v , we must account for ad- 
ditional uncertainty: viz 
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Tt = 3072.3 ±0.9 ± 1.1, (4) 

where the first error is the statistical error of the fit (as illustrated in Fig. Q), 
and the second is an error related to the systematic difference between the two 
calculations of 5c by Towner, Hardy and HarveyB and by Ormand and BrownEl] 
that we have combined in reaching this result. (For a more complete discussion 
of how we treat these theoretical uncertainties, see referenced.) We now add 
the two errors linearly to obtain the value we use in subsequent analysis: 



.Fi = 3072.3 ± 2.0. (5) 

The value of V u d is obtained by relating the vector constant, G v , deter- 
mined from this Tt value, to the weak coupling constant from muon decay, 
G F /(hcf = (1.16639 ±0.00001) x 10~ 5 GeV~ 2 , according to: 

2 K 

Vud = 2GU1 + ADW (6) 

The result obtained is 



| Kdh 0.9740 ±0.0005, [Nuclear] (7) 

where the nucleus-independent radiative correction has been set at 

= (2.40 ±0.08)%. (8) 

Note this value differs slightly £tjut within errors) from an earlier value 
because of the decision by Sirlinli3 to centre the cut-off parameter m A , where 
( TO ai/2) < m A < 2m ai , exactly at the ai-meson mass when evaluating the 
axial contribution to the radiative-correction loop graph. 

From the value of V u d given in Eq. (]?]), the unitarity sum, Eq. (^|), becomes 

= 0.9968 ± 0.0014, [Nuclear] (9) 

i 

which fails to meet unity by 2.2 standard deviations. In connection with this 
result, we note the following two points: 

(a) The error bar associated with | V u d | in Eq. (^) is not predominantly 
experimental in origin. In fact, if experiment were the sole contributor, the 
uncertainty would be only ±0.0001. The largest contributions to the | V u d | 
error bar come from (±0.0004) and S c (±0.0003). 

(b) The unitarity result in Eq. (|^) depends on the values of nuclear structure- 
dependent corrections. In Sec. S we will examine whether the failure to meet 
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unitarity can be repaired by reasonable adjustments to these corrections. Our 
conclusion is largely negative. Other speculative possibilities are presented in 
Sec. | 



2.2 Neutron decay 

On the one hand, free neutron decay has an advantage over nuclear decays 
since there are no nuclear-structure dependent corrections to be calculated. 
On the other hand, it has the disadvantage that it is not purely vector-like 
but has a mix of vector and axial-vector contributions. Thus, in addition to 
a lifetime measurement, a correlation experiment is also required to separate 
the vector and axial-vector pieces. Both types of experiment present serious 
experimental challenges. The value of V u d is determined from the expression 

V 2 = K/k2 

ud G?(l + Av)(l + 3A a )/(l + « fl )T n ' 1 ' 

where A is the ratio of axial- vector and vector effective coupling constants, 
A = G'JG' V , with G'l = + A*) and G'^ = G%{1 + AJQ. Here A* and 
A^ are the nucleus-independent radiative corrections. With the experimental 
measurement of A being a determination of G' A /G' V , the actual value of A* is 
not required for the evaluation of V^ d . In Eq. (|l0|), / is the statistical rate 
function and 5r, the nucleus^dependent radiative correction evaluated for the 
case of a neutron. WilkinsonEa has,evaluated the product f(l + Sr). His value 
was revised by Towner and HardyEJ who incorporated the current best Q- value 
to obtain /(l + S R ) = 1.71489 ± 0.00002. Lastly, r„ is the mean lifetime for 
neutron decay. 

A survey of world data on neutron decay appears in Table ||. The life- 
time measurements, when considered as a single body of data, are statistically 
consistent. However the measurements of A are not (x 2 jv = 2.0) and, as a 
consequence, the uncertainty quoted in the table for the overall average value 
of A has been scaled by a factor of 1.7. Inserting the average values for r„ and 
A into Eq. (|l(]), we determine the value of V u d to be 

| Kd |= 0.9759 ±0.0021, [Neutron] (11) 

and the unitarity sum to be 

Y^V^ = 1.0007 ±0.0042, [Neutron] (12) 
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Table 2: Experimental results" for neutron decay 



Method 



Measured values 



Average 



j(s) n beam 



n trap 



918 ± 14 




891 ± 90-, 




876 ± 21 




878 ± 3lllf 




889.2 ±i 




877 ± loE> 




903 ± 13 






887.6 ±3.0E 


888.4 ±3.3 


22 


882.6 ±2.7E 


885.4 ± 1.0 


24 



Overall Average 

/3-asym. -1.254 ± 0.015^ 
-1.262 ±0.005^ 
-1.266 ±0.004 22 

e-Fcorr. -1.259 ±0.017 

Overall Average 



33 



-1.257 ±0.012 ft [ 
-1.2594 ±0.00»8J 291 
-1.274 ±0.003 1? ^ 



891.2 ±4.8 



885.5 ±0.9 
885.8 ±0.9 



-1.2665 ±0.0031 c 
-1.259 ± 0.017 
-1.2664 ±0.0031 



a Following the practice used for superallowed decay, we retain only those measurements 
with uncertainties that are within a factor of ten of the most precise measurement for each 
quantity. All such measurements, of which we are aware, that have not been withdrawn by 
their authors are listed. 1— 1 

b Corrected for weak magnetism and recoil following ref.cj . 1—1 
c Not included is a preliminary value of —1.2735 ± .0021 reported at this conference E3 . 



a value that agrees with unitarity and with the nuclear result, Eq. (^), which 
is a factor of three more precise. In connection with this result, we note the 
following two points: 

(a) For neutron decay, the error bar associated with | V u d | in Eq. (|ll|) is 
some four times larger than the error bar obtained from nuclear decays, Eq. 

; however, in contrast with the latter case, it is predominantly experimental 
in origin. The largest theoretical contribution to the | V u d | error bar comes (at 
the level of ±0.0004) from A^, a correction that is common to both nuclear 
and neutron decays. 

(b) Currently, the theoretical uncertainty on dominates the nuclear 
result for | V u d |. As experimental results for the neutron improve, A^ will 
eventually dominate the neutron result too. Therefore, so long as A^ remains 
at its current level of uncertainty, the neutron results will never be able to test 
unitarity with significantly better precision than the nuclear decays do now, 
in spite of their independence of 8c- They will, of course, be able to provide 
an important test of whether there are some systematic problems with the 
nuclear-dependent corrections, which are not now anticipated in the theoretical 
uncertainty quoted in Eq. (pLL 

At the conference, Reichl23 presented a new result for the beta asymmetry 
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obtained by the PERKEO II collaboration. This result is more precise than 
any of its predecessors and leads to a value of A = —1.2735 ± 0.0021. If we 
combine this single value - i.e., not averaged with its predecessors - together 
with the current world average for the neutron lifetime, Eq. (foh, yields 



\V ud \= 0.9714 ± 0.0015, [Neutron PERKEO II] (13) 

and a unitarity sum of 

Vut = 0-9919 ± 0.0030. [Neutron PERKEO II] (14) 

i 

The error bars here are smaller than they were for the world-average neutron 
results, although | V u d | is still dominated by the uncertainty in the beta asym- 
metry measurement. The unitarity sum itself, though, is tantalizingly similar 
to the nuclear result in that it is less than unity by several standard deviations. 

2.3 Pion beta decay 

Like neutron decay, pion beta decay has an advantage over nuclear decays 
in that there are no nuclear structure-dependent corrections to be made. It 
also has the same advantage as the nuclear decays in being a purely vector 
transition, in its case 0~ — > 0~, so no separation of vector and axial- vector 
components is required. Its major disadvantage, however, is that pion beta 
decay, ir + — > ir e + v e , is a very weak branch, of the order of 10~ 8 . This results 
in severe experimental limitations. For the pion beta decay, the value of V u d 
is determined from the expression 

ud Gl{l + Al)hhf{l + 8 R )T m /B& [ ° ! 

where / is the approximate statistical rate function 

/-if^V (16) 



30 \m , 

with A being the pion mass difference (A = m w + — m T o) and m e , the electron 
mass. The mass difference is known with high precision from the work of 
Crawford et aM to be A = (4.5936 ± 0.0005) MeVp-The factors fx and f 2 are 
corrections to /, and are easily calculable functionsE3 of A/m„.+ with values of 
fi = 0.941039 and f 2 — 0.951439. The nucleus-dependent radiative correction 
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evaluated for the case of the pion is Sr = (1.05 ± 0.15)%. Finally, r m and BR 
are the pion mean lifetime and branching ratio, respectively. 

Two precise lifetime measurements E3'E3 were published in 1995 and the 
PDG98 average is 

r m = (2.6033 ± 0.0005) x 10~ 8 s. (17) 
The branching ratio is from McFarlane et a/S: 

BR = (1.025 ±0.034) x 10~ 8 . (18) 
Inserting Eqs. (|l7|) and @ into Eq. (jl|), we determine the value of V u d to be 

| Kd|= 0.9670 ±0.0161, [Pion] (19) 

and the unitarity sum 

= 0.9833 ± 0.0311, [Pion] (20) 

i 

satisfying the unitarity condition but with comparatively large uncertainty. 
The error on LV u d I is entirely due to the uncertainty in the pion branching 
ratio. Deutsclio reported to this conference that there is a proposal submitted 
to PSI for an experiment to improve this branching ratio by a factor of eight. 
If successful, the error on | V u d | would be reduced to ±0.0022, comparable 
to the current limit from neutron decay. Of course, ultimately it too will be 
limited by the theoretical uncertainty on 



3 Theoretical corrections in nuclear decays 



In Sec. 2.1 we noted that the value of V u d determined from nuclear decays 
- the most precise result available - resulted in the unitarity test among the 
elements of the first row of the CKM matrix not being satisfied by two stan- 
dard deviations. Here, we discuss the theoretical corrections involved in the 
determination, and assess whether the failure to meet unitarity can be removed 
by reasonable adjustments in these calculations. To restore unitarity, the cal- 
culated radiative corrections for all nine nuclear transitions would have to be 
shifted downwards by 0.3% (i.e. as much as one-quarter of their current value), 
or the calculated Coulomb correction shifted upwards by 0.3% (over one-half 
their value), or some combination of the two. 
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Table 3: Calculated nucleus-dependent radiative correction, Sr, in percent units, and the 
component contributions as identified in Eq. 







£ 5 2 


^3 






Sr 










unquenched 


quenched 


quenched 


io c 


1.47 


0.18 


0.01(1) 


-0.39(5) 


-0.36(4) 


1.30(4) 


14Q 


1.29 


0.23 


0.01(1) 


-0.27(7) 


-0.26(5) 


1.26(5) 


26m Al 


1.11 


0.33 


0.02(2) 


0.06(1) 


-0.01(1) 


1.45(2) 


34 C1 


1.00 


0.39 


0.03(3) 


-0.04(1) 


-0.09(1) 


1.33(3) 


38m |£ 


0.96 


0.41 


0.04(4) 


-0.02(2) 


-0.09(2) 


1.33(4) 


42 Sc 


0.94 


0.45 


0.05(4) 


0.12(2) 


0.03(2) 


1.47(5) 


46y 


0.90 


0.47 


0.06(6) 


0.04(1) 


-0.03(1) 


1.40(6) 


50 Mn 


0.87 


0.49 


0.07(7) 


0.04(1) 


-0.03(1) 


1.40(7) 


54 Co 


0.84 


0.51 


0.07(7) 


0.05(1) 


-0.03(1) 


1.40(7) 



3.1 Radiative corrections 
As mentioned in Sec 



2.1 



the radiative correction is conveniently divided into 
terms that are nucleus-dependent, 5r, and terms that are not, A^. These are 
written 



a 
a 



[g(E m )+6 2 + S 3 + 2C NS ] 



K = 7T [41n(m z /m p ) + ln(m p /m A ) + 2C E 

Z7T 



(21) 



where the ellipses represent further small terms of order 0.1%. In these equa- 
tions, E m is the maximum electron energy in beta decay, m z the Z-boson 
mass, m A the ai-meson mass, and 62 and ^3 the order- ifa 2 and -Z 2 a 3 contri- 
butions respectively. The function_n(i? e , E m ), which is a function of electron 
energy, was first defined by Sirlin U2I as part of the order-a universal photonic 
contribution arising from the weak vector current; it is here averaged over the 
electron spectrum to give g{E m ). Finally, the terms Ceom and Cns come 
from the order-a axial-vector photonic contributions: the former accounts for 
single- nucleon contributions, while the latter covers two-nucleon contributions 
and is consequently dependent on nuclear structure. 

Calculated values for all four components of Sr are given in Table ||. There 
have been two independent calculations e!Ie3 of both 82 and £3; they are com- 
pletely consistent with one another if proper account is taken of finite-size 
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effects in the nuc ear charge distribution. The._values listed in Table |3| are 
our recalculations □ using the formulas of SirlinEU but incorporating a Fermi 
charge-density distribution for the nucleus. Note that we have followed Sirlin 
in assigning an uncertainty equal to {a/2-K)5^, as an estimate of the error made 
in stopping the calculation at that order. Also appearing in the table are two 
values of (cc/i^Cns for each decay, according to whether the weak axial and 
electromagnetic couplings at the nucleoli take their free-nucleon values CJ (un- 
qucnchcd) or medium-modified valuesCj (quenched). We adopt the quenched 
values in evaluating 5r. 

In assessing the changes in Sr that would be required in order to restore 
unitarity, it is helpful to rewrite Eq. ( pl| ) in terms of the typical values taken 
by its components: viz 

S R ~ 1.00 + 0.40 + 0.05 + (a/ir)C NS %, (22) 

where (o./it)Cns is of order —0.3% for T z = —1 beta emitters, 10 C and 14 0, 
and of order five times smaller for the T z = emitters, ranging from —0.09% 
to +0.03%. Thus, for T z = emitters S R ~ 1.4%. If the failure to obtain 
unitarity in the CKM matrix with V u d from nuclear beta decay is due to the 
value of 5r, then 5r must be reduced to 1.1%. This is not likely. The leading 
term, 1.00%, involves standard QED and is well verified-. The order-Za 2 term, 
0.40%, while less secure has been calculated twice EJc3 independently with 
results in accord. 

Taking a similar approach for the nucleus-independent radiative correction, 
we write 



A^ = 2.12 - 0.03 + 0.20 + 0.1% ~ 2.4%, (23) 

of which the first term, the leading logarithm, is unambiguous. Again, to 
achieve unitarity of the CKM matrix, would have to be reduced to 2.1%: 
i.e. all terms other than the leading logarithm must sum to zero. This also 
seems unlikely. 



3.2 Coulomb corrections 

Because the leading terms in the radiative corrections are so well founded, at- 
tention has focussed more on the Coulomb correction. Although smaller than 
the radiative correction, the Coulomb correction is clearly sensitive to nuclear- 
structure issues. It comes about because Coulomb and charge-dependent nu- 
clear forces destroy isospin symmetry between the initial and final states in 
superallowed beta-decay. The consequences are twofold: there are different 
degrees of configuration mixing in the two states, and, because their binding 
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Table 4: Calculated Coulomb correction, 8q, in percent units. 



THH a OB b SVS C NBO rf Adopted 6 

Value 



10 C 0.18 0.15 0.00 0.12 0.16(3 

14 0.28 0.15 0.29 0.22(3 

™A1 0.33 0.30 0.27 0.31(3 

;4 C1 0.64 0.57 0.33 0.61(3 

m K 0.64 0.59 0.33 0.62(3 

l2 Sc 0.40 0.42 0.44 0.41(3 

46 V 0.45 0.38 0.41(3 

Mn 0.47 0.35 0.41(3 

4 Co 0.61 0.44 0.49 0.52(3 



flm firim rnfs. BEl : frnn from ref. R . 
b Ref. El ; c Ref. H ; d Ref. 1*3 . 

e Average of OB and THH; assigned uncertainties reflect the 
relative scatter between these calculations. 



energies are not identical, their radial wave functions differ. Thus, we accom- 
modate both effects by writing 5c = Sci + &ci- Constraints can be placed on 
the calculation of 8ci by insisting that it reproduce the measured coefficients 
of the isobaric mass multiplet equation. Constraints are placed on 5ci by in- 
sisting that the asymptotic forms of the proton and neutron radial functions 
match the known separation energies. 

The results of several calculations for 5c are shown in Table ||. The values 
in the first column_arc those calculated by the methods developed by Towner, 
Hardy and Harveyu and refined in more recent publications! 4 !]'! 4 !!. They result 
from shell- model calculations to determine 5ci, and full-parentage expansions 
in terms of Woods-Saxon radial wave functions to obtain 5c2- Ormand and 
Brown, whose values E3 for 5c appear in column 2, also employed the shell 
model for calculating 5ci but derived 5c2 from a self-consistent Hartree-Fock 
calculation. Both of these independent calculations for 5c reproduce the mea- 
sured coefficients of the relevant isobaric multiplet mass equation, the known 
proton and neutron separation energies^ and the measured ft- values of the 
weak non-analogue + — > + transitions! 4 !! where they are known. In our anal- 
ysis in Sec. 2.1, we have used the average of these two sets of 5c values: our 
adopted values appear in the last column of the table. 

Two more recent calculations provide a valuable check that these 5c values 
are not suffering from severe systematic effects. Sagawa, van Giai and SuzukieS 
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have added RPA correlations to a Hartree-Fock calculation that incorporates 
charge-symmetry and charge-independence breaking forces in the mean-field 
potential to take account of isospin impurity in the core; the correlations, 
in essence, introduce a coupling to the isovector monopole giant resonance. 
The calculation is not constrained, however, to reproduce known separation 
energies. Finally, a large shell-model calculation has been mounted for the 
A = 10 case by Navratil, Barrett and Ormandca. Both of these two new works 
have produced values of 5c very similar to, but actually smaller than those 
used in our analysis, i.e. worsening rather than helping the unitarity problem. 

The typical value of 5c is of order 0.4%. If the unitarity problem is to be 
solved by improvements in 5c, then 5c has to be raised to around 0.7%. There 
is no evidence whatsoever for such a shift from recent works. 

4 Speculative suggestions to resolve unitarity problem 

Since it is concluded in Sec. [3] that reasonable adjustments to the theoretical 
corrections, Sr, AJ( and 5c, are unlikely to resolve the unitarity problem posed 
by the nuclear result for V u d , we turn in this section to some more speculative 
alternatives. We discuss four suggestions, two that do not require extensions 
to the Standard Model and two that do. All but one require the introduction 
of important new physics to explain the apparent discrepancy. 

4-1 Saito- Thomas correction 

One suggestion that would resolve much of the unitarity problem was suggested 
by Thomas at the last WEIN conference and subsequently published as a let- 
ter l3. It is based on a quark- meson coupling model, in which nuclear matter 
consists of non-overlapping nucleon (MIT) bags bound by the self-consistent 
exchange of a and to mesons in the mean-field approximation. The model is ex- 
tended to include an isovector-vector meson (/?) and an isovector-scalar meson 
(5). The coupled-equations to be solved are very similar to those of Quantum 
Hadrodynamics Ej, but involve boundary conditions at the bag radius. As a 
consequence of a coupling to meson fields the quark mass in a medium becomes 
an effective one, 

m* =m,- (V„ ± \V S ), i = u,d (24) 

with the upper sign for the up quark. Here V a and Vg are strengths of a- 
meson and <5-meson mean fields. At the quark level, the conserved vector 
current (CVC) hypothesis is broken if the up and down quarks have different 
masses, but for a free nucleon this breaking is second order in (m<2 — m u ). In a 
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nuclear medium, however, this breaking may become first order in (mi — to*), 
but a critical constraint of the calculation is to show it reverts to (m^ — m u ) 2 
at zero density. 

Saito-Thomas (ST) solve the coupled equations in an infinite nuclear- 
matter approximation, the solutions being functions of the matter density, 
Pb- Let ifii/j be the wavefunction of quark i bound in nucleon j in nuclear 
matter; then, the following overlap integral can be defined between quark i 
bound in a proton and quark i' bound in a neutron: 

Iw{pb)= [ dV^.A'/n- (25) 

JBag 

Calculations indicate that Iu'(pb) varies linearly in pb, being unity at zero 
density. What is required for beta decay is the product of three overlap inte- 
grals 

\I ud \ 2 x\I uu \ 2 x\I dd \ 2 ^l-8^ (26) 

and this product also is approximately linear in the density. The results of the 
calculations E°l are 

^ uark = b x , (27) 

where po is the saturation density for equilibrium nuclear matter, and b is in 
the range (0.15 — 0.20)% for bag radii in the range (0.6 — 1.0) fm. Thus <5p Ualk , 
at densities of po/2, lies in the range 0.075% to 0.10% while, at po, it lies 
between 0.15% and 0.20%. The corresponding value for V 2 d would then be 
increased by these amounts. 

It is clearly premature to apply such a correction at this stage of its de- 
velopment, although the Particle Data Group has done so in its most recent 
surveyH. The authors themselves admitLJ that their results are merely quali- 
tative, having been derived from a model that deals only with nuclear matter. 
Quantitative results must await extensions of the formalism so that it can be 
applied in finite nuclei. There is also a question of whether simply adding the 
Saito-Thomas correction of <5ji, uark to the Sq already computed in a nuclcons- 
only calculation is the correct approach at all, since it may lead to double 
counting. In the nucleons-only case, certain parameters are adjusted to re- 
produce Coulomb observables such as the b- and c-coefficients of the isobaric 
mass multiplet equation, and proton and neutron separation energies. Thus, 
the parameters become effective ones- which in some unquantifiable way, actu- 
ally contain quark effects. ThomasE3 disagrees, however, claiming that there 
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is no double counting because the nuclear shell model is derivable from the 
quark-meson coupling model. 

However imprecise 5^ uark may be, it must be admitted that it would act to 
reduce the unitarity problem, supplying possibly up to 0.2% of the 0.3% dis- 
crepancy. This is where the free-neutron decay will provide a decisive answer, 
since it requires no Saito-Thomas correction. The essence of the model is the 
coupling of quarks to cr-mesons, and these scalar fields are only present in a 
medium. 

Finally, if this correction turns out to be warranted in the nuclear case, it 
would be an extremely interesting result in its own right even though it would 
obviate the need for an extension to the Standard Model. Such a result would 
constitute the first case in which quark degrees of freedom were unequivocally 
required to explain a nuclear-structure observable. 

4-2 Ad hoc two-parameter fit 

Wilkinson @ has suggested that data like those displayed in Fig. [I] might be 
better fitted by a two-parameter function 

Tt = Tt(0) [1 + aZ] , (28) 

where the term proportional to Z represents a further correction of unknown 
origin. A fit of such a function to the data yields 



a = (0.77 ±0.55) x 10~ 4 
Tt(0) = (3068.3 ± 3.2) s (29) 

with x 2 /v — 0.86. In this case, the error is dominated by the statistical error of 
the fit, with only a small contribution due to the systematic difference between 
the two calculations of 6c (see the discussion following Eq. (Q)). If this value of 
Tt(0) is used in Eq. @, the unitarity sum, (Eq. (§)), becomes 0.9981 ±0.0016, 
which is acceptably close to one. 

Although this result for the unitarity test is less provocative than the one 
presented in Eq. (^|), we emphasize that there is no physical justification what- 
soever for incorporating any Z-dependent corrections beyond those already 
accounted for in 5r and 8c- Furthermore, there is no statistically significant in- 
dication from the Tt- value data that one is required. Consequently, until such 
time as new Tt-vahie data can demonstrate a clear residual Z-dependence, 
this suggested explanation for the nuclear unitarity result must be regarded as 
unsatisfactory. 
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4-3 Right-hand currents 



Because experimental evidence at low energies favors maximal parity violation 
in weak interactions, that condition has been built into the Standard Model. 
However, one class of possible extensions to the Standard Model would re- 
store parity symmetry at higher energies through the introduction of additional 
heavy charged gauge bosons that are predominantly right-handed in character. 
Under certain specific conditions, such models could explain the nuclear results 



described in Sec. 2.1 



Fox_example, an extension._known as the manifest left-right symmetric 
model til leads to a revised formEJ for Eq. (|^) which includes a mixing angle £: 

tte-'O- ManW (30) 

If, under these conditions, we require that V 2 d must satisfy unitarity, viz 

V* d = 1 - Vl ~ V 2 b , (31) 

then, with the Tt value taken from Eq. (^), we can derive a value for the 
mixing angle of £ = 0.0015 ± 0.0007. Within the context of the manifest left- 
right symmetric model, this is a fully satisfactory result; not surprisingly, it is 
two standard deviations away from the Standard Model's pure V — A value. 



4-4 Scalar interaction 

If the Standard Model of pure V—A weak interactions is extended instead to in- 
clude scalar and tensor interactions, then the expression for the beta spectrum 
shape would include an additional term coming from the scalar interaction, 
which is inversely proportional to the electron energy. Since the Tt values 
include an integral over the spectrum shape, the presence of a non-zero scalar 
interaction would be reflected by a correction to the Tt values that is inversely 
proportional to the decay energy of the parent nucleus. To test this possibility, 
we fit the data in Fig. [l] by the function 

Tt = k[l + b Fl (W- 1 )} , (32) 

where k and bp are parameters determined from the fit. The latter, bp, is 
known as the Fierz interference term; a non-zero value for this term would 
signal the presence of a scalar interaction. In Eq. (j32|), -f 2 = (1 — (aZ) 2 ), a 
is the fine-structure constant, Z the atomic number of the daughter nucleus, 
and is the value of l/W averaged over the electron spectrum, where W 
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is the electron energy in electron rest-mass units. The value obtained for the 
Fierz interference term from the data in Fig. |l| is 

b F = -0.034 ±0.0026. (33) 

A negative sign is expected @ for a positron emitter. The 90% confidence level 
upper limit is 

\b F \< 0.0077. (34) 

To proceed to V u d, we need a value of J-t extrapolated to the Z = limit. To 
this end, we fit the nine values of 7(W _1 ) by a second-order polynomial in Z 
and use this polynomial to provide the value of 7(M /_1 ) at Z = 0. With k and 
b F taken from the fit, Eq. ((33|), we obtain 

Tt(0) = 3066.4 ±3.1 

\Vud\ = 0.9749 ± 0.0006 [Nuclear + bF] 

V*i = 0.9986 ± 0.0016. [Nuclear + bF] (35) 

i 

This result is now in accord with unitarity but, of course, requires the presence 
of a non-zero scalar coupling to achieve that goal. 

We close this section by noting that all four suggested explanations we 
have presented for the non-unitarity result in Eq. (^]) are entirely speculative 
and should be considered with caution. 

5 Future experimental prospects 

It is evident from the foregoing discussion that the current world data on 
V u d are tantalizingly close to producing a definitive result on unitarity. The 
nuclear measurements have achieved the highest experimental precision but 
they are now constrained by theoretical uncertainties. The neutron and pion 
measurements are, as yet, experimentally less precise than the nuclear ones, but 
they are free from one of the more important sources of theoretical uncertainty, 
8c- All three classes of measurements are now being extended and improved 
at a number of laboratories, and there are good prospects for considerably 
reduced error iars within a few years. In combination with the re- visitation of 
the K e % decayta, these results should settle the uncertainty over 5c and determie 
whether the deviation from unitarity, apparent from the nuclear result, is real 
or not. This is ample argument to justify considerable experimental activity. 
At the same time, it should not be forgotten that the full impact of these 
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experimental advances will be diluted until there are theoretical improvements 
in the correction - the dominant source of theoretical uncertainty in all 
cases. Ultimately, to make significant improvements in the unitarity test, there 
will have to be advances in both theory and experiment. 

Of the three classes of experiment, that focusing on pion decay is currently 
farthest from the precision required for a meaningful unitarity test. We noted 
in Sec. 2.3 that a considerable improvement is anticipated in the foreseeable 
future, but the result is still not likely to reach a precision higher than that 
already achieved for neutron decay. There is considerable optimism, however, 
that the neutron measurements themselves can be improved as new experi- 
ments with ultra-cold neutrons come to fruition. The outcome can be expected 
seriously to challenge the nuclear experiments for experimental precision, but 
will take at least a few more years to do so. 

As to the nuclear experiments, the nine superallowed transitions whose 
ft values are known to within a fraction of a percent have been the subject 
of intense scrutiny for at least the past three decades. All except 10 C have 
the special advantage that the superallowed branch from each is by far the 
dominant transition in its decay (> 99%). This means that the branching 
ratio for the superallowed transition can be determined to high precision from 
relatively imprecise measurements of the other weak transitions, which can 
simply be subtracted from 100%. Given the quantity of careful measurements 
already published, are there reasonable prospects for significant improvements 
in these decay measurements in the near future? Given the uncertainty in 
the theoretical corrections, which experiments can shed the most light on the 
efficacy of these corrections? 

If we begin by accepting that it is valuable for experiment to be at least 
a factor of two more precise than theory, then an examination of the world 
data shows that the Q-values for 10 C, 14 O, 26m Al and 46 V, the half-lives of 
10 C, 34 C1 and 38m K, and the branching ratio for 10 C can all bear improve- 
ment. Such improvements will soon be feasible. The Q-values will reach the 
required level (and more) as mass measurements with new on-line Penning 
traps become possible; half-lives will likely yield to measurements with higher 
statistics as high-intensity beams of separated isotopes are developed for the 
new radioactive-beam facilities; and, finally, an improved branching-ratio mea- 
suremeni-pn 10 C has already been made with Gammasphere and simply awaits 
analysis E3. 

Qualitative improvements will also come as we increase the number of 
superallowed emitters accessible to precision studies. The greatest attention 
recently has been paid to the T z — emitters with A > 62, since these nu- 
clei are expected to be produced at new radioactive-beam facilities, and their 
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calculated Coulomb corrections, 5c, are predicted to be larg eBBEI. They 
could then provide a valuable test of the accuracy of 5c calculations. It is 
likely, though, that the required precision will not be attainable for some time 
to come. The decays of these nuclei will be of higher energy and each will 
therefore involve several allowed transitions of significant intensity in addition 
to the superallowed transition. Branching-ratio measurements will thus be 
very demanding, particularly with the limited intensities likely to be available 
initially for these rather exotic nuclei. Lifetime measurements will be similarly 
constrained by statistics. 

More accessible in the short term will be the T z = — 1 superallowed emit- 
ters with 18 < A<~ 38. There is good reason to explore them. For example, the 
calculated value □ of 5c for 30 S decay, though smaller than the Sc's expected 
for the heavier nuclei, is actually 1.2% - about a factor of two larger than for 
any other case currently known - while 22 Mg has a very low value of 0.35%. 
If such large differences are confirmed by the measured /t-values, then it will 
do much to increase our confidence in the calculated Coulomb corrections. To 
be sure, these decays will provide a challenge, particularly in the measurement 
of their branching ratios, but the required precision should be achievable with 
isotope-separated beams that are currently available. In fact, such experiments 
are already in their early stages at the Texas A&M cyclotron. 

6 Conclusions 

The current world data on superallowed + — » + beta decays lead to a 
self-consistent set of .Ft- values that agree with CVC but differ provocatively, 
though not yet definitively, from the expectation of CKM unitarity. There are 
no evident defects in the calculated radiative and Coulomb corrections that 
could remove the problem, so, if any progress is to be made in firmly establish- 
ing (or eliminating) the discrepancy with unitarity, additional experiments are 
required. We have indicated what some relevant nuclear experiments might 
be. 

In the past decade there have been significant improvements in the mea- 
surements of the neutron lifetime and beta asymmetry, and further improve- 
ments are promised in the near future. It is likely that these studies in neutron 
decay will soon approach the results from nuclear superallowed decays in pre- 
cision; and when there, they will have the advantage in that their results are 
not dependent on nuclear-structure corrections. Intriguingly, the most recent 
result E3 from PERKEO II is yielding a low value of V u d and a failure by two 
standard deviations of CKM unitarity: a result the nuclear decays have had 
for several years. We have offered some speculative suggestions as to what the 
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nuclear failure could be due to. If the neutron results exhibit the same failure, 
though, only the suggestions involving extensions to the Standard Model could 
possibly apply. 

Clearly, there is strong motivation to pursue experiments on both the neu- 
tron and nuclear front, since, if firmly established, a discrepancy with unitarity 
would indicate the need for important new physics. 

Acknowledgments 

We would like to thank Tony Thomas for comments on the Saito-Thomas cor- 
rection. 1ST acknowledges the hospitality of the Theoretical Physics Division 
of the Los Alamos National Laboratory, where he was when this report was 
written. The work of JCH was supported by the U.S. Department of En- 
ergy under Grant number DE-FG05-93ER40773 and by the Robert A. Welch 
Foundation. 

References 

1. N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963) 

2. M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973) 

3. C. Caso et al, Eur. Phys. J. C 3, 1 (1998) 

4. W.J. Marciano, these proceedings 

5. J.C. Hardy, I.S. Towner, V.T. Koslowsky, E. Hagberg and H. Schmeing, 
Nucl. Phys. A 509, 429 (1990) 

6. I.S. Towner, E. Hagberg, G. Savard, J.C. Hardy, V. Koslowsky, A. 
Galindo-Uribarri and D.C. Radford, Proceedings of the IV Int. Symp. on 
Weak and Electromagnetic Interactions in Nuclei, Osaka, Japan, June 
1995, eds. H. Ejiri, T. Kishimoto and T. Sato (World Scientific, Singa- 
pore, 1995) pp. 313-318 

7. B.K. Fujikawa et al, APS e-print 1998may29_004 

8. P.D. Harty, N.S. Bowden, P.H. Barker and PA. Amundsen, Phys. Rev. 
C 58, 821 (1998) 

9. I.S. Towner, J.C. Hardy and M. Harvey, Nucl. Phys. A 284, 269 (1977) 

10. W.E. Ormand and BA. Brown, Phys. Rev. C 52, 2455 (1995) 

11. W.J. Marciano and A. Sirlin, Phys. Rev. Lett. 56, 22 (1986) 

12. A. Sirlin, in Precision Tests of the Standard Electroweak Model, ed. P. 
Langacker (World-Scientific, Singapore, 1994) 

13. C.J. Christenscn, A. Neilsen, A. Bahnsen, W.K. Brown and B.M. Rustad, 
Phys. Rev. D 5, 1628 (1972) 

14. P.E. Spivak, Sov. Phys. JETP 67, 1735 (1988) 



20 



15. J. Last, M. Arnold, J. Dohner, D. Dubbers and S.J. Freedman, Phys. 
Rev. Lett. 60, 995 (1988) 

16. R. Kossakowski, P.Grivot, P. Liaud, K. Schreckenbach and G. Azuelos, 
Nucl. Phys. A 503, 473 (1989) 

17. J. Byrne, P.G. Dawber, C.G. Habeck, S.J. Smidt, J.A. Spain and A.P. 
Williams, Europhys. Lett. 33, 187 (1996) 

18. Yu.Yu. Kosvintsev, V.I. Morozov and G.I. Terekhov, JETP Lett. 44, 
571 (1986) 

19. V.I. Morozov, Nucl. Instrum. Methods A 284, 108 (1989) 

20. W. Paul, F. Anton, L. Paul, S. Paul and W. Mampe, Z. Phys. A 45, 25 
(1989) 

21. W. Mampe, P.Ageron, J.C. Bates, J.M. Pendlebury and A. Steyerl, Phys. 
Rev. Lett. 63, 593 (1989) and Nucl. Instrum. Methods A 284, 111 (1989) 

22. V.V. Nesvizhevskii, A.P. Scrcbrov, R.R. Tal'dacv, A.G. Kharitonov, V.P. 
Alnmenkov, A.V. Strclkov and V.N. Shvctsov, Sov. Phys. JETP 75, 
405 (1992) 

23. W. Mampe, L.N. Bondarenko, V.I. Morozov, Yu. N. Panin and A.I. 
Fomin, JETP Lett. 57, 82 (1993) 

24. P. Geltenbort, WEIN98; S. Arzumanov, L. Bondarenko, S. Chernyavsky, 
W. Drexel, A. Fomin, P. Geltenbort, F. Gonncnwein, V. Morozov, Yu. 
Panin, J. Pendlebury and K. Schreckenbach, in V Int. Seminar on Inter- 
action of Neutrons with Nuclei, Dubna 1997 pp. 53-66 

25. V.E. Krohn and G. R. Ringo, Phys. Lett. B 55, 175 (1975) 

26. D.H. Wilkinson, Nucl. Phys. A 377, 474 (1982) 

27. B.G. Erozolimskii, A.I. Frank, Yu. A. Mostovoi, S.S. Arzumanov and 
L.R. Voitzck, Sov. J. Nucl. Phys. 30, 356 (1979) 

28. P. Bopp, D. Dubbers, L. Hornig, E. Klcmt, J. Last, H. Schiitze, S.J. 
Freedman and O. Scharpf, Phys. Rev. Lett. 56, 919 (1986) 

29. B. Yerozolimsky, I. Kuznetsov, Yu. Mostovoy, I. Stepanenko, Phys. Lett. 
B 412, 240 (1997) 

30. P. Liaud, K. Schreckenbach, R. Kossakowski, H. Nastoll, A. Bussiere, 
J.P. Guillaud and L. Beck, Nucl. Phys. A 612, 53 (1997) 

31. H. Abele, S. Baessler, D. Dubbers, J. Last, U. Mayerhofer, C. Metz, 
T.M. Miillcr, V Nesvizhevsky, C. Raven, O. Scharpf and O. Zimmer, 
Phys. Lett. B 407, 212 (1997) 

32. J. Reich, WEIN98; J. Reich, H. Abele, M. Astruc Hoffmann, S. Baessler, 
P. v. Biilow, D. Dubbers, V. Nesvizhevsky, U. Peschke and O. Zimmer, 
in V Int. Seminar on Interaction of Neutrons with Nuclei, Dubna 1997 

33. C. Stratowa, R. Dobrozemsky and P. Weinzierl, Phys. Rev. D 18, 3970 
(1978) 



21 



34. I.S. Towner and J.C. Hardy, in Symmetries and Fundamental Interac- 
tions in Nuclei, eds. E.M. Henley and W.C. Haxton (World-Scientific, 
Singapore, 1995) pp.183-249 

35. J.F. Crawford, M. Damn, R. Frosch, B. Jost, P.-R. Kettle, R.M. Marshall, 
B.K. Wright and K.O.H. Ziock, Phys. Rev. D 43, 46 (1991) 

36. W.K. McFarlane et al, Phys. Rev. D 32, 547 (1985) 

37. V.P. Koptev et al, JETP Lett. 61, 877 (1995) 

38. T. Numao, J. A. Macdonald, G.M. Marshall, A. Olin and M.C. Fujiwara, 
Phys. Rev. D 52, 4855 (1995) 

39. J. Deutsch, these proceedings 

40. A. Sirlin, Phys. Rev. 164, 1767 (1967) 

41. A. Sirlin, Phys. Rev. D 35, 3423 (1987); A. Sirlin and R. Zucchini, Phys. 
Rev. Lett. 57, 1994 (1986) 

42. W. Jaus and G. Raschc, Phys. Rev. D 35, 3420 (1987) 

43. I.S. Towner, Nucl. Phys. A 540, 478 (1992) 

44. I.S. Towner, Phys. Lett. B 333, 13 (1994) 

45. W.E. Ormand and B.A. Brown, Phys. Rev. Lett. 62, 866 (1989); Nucl. 
Phys. A 440, 274 (1985) 

46. I.S. Towner, in Symmetry Violations in Subatomic Physics, eds. B. Castel 
and P.J. O'Donnel (World-Scientific, Singapore, 1989) p. 211 

47. E. Hagberg, V.T. Koslowsky, J.C. Hardy, I.S. Towner, J.C Hykawy, G. 
Savard and T. Shinozuka, Phys. Rev. Lett. 73, 396 (1994) 

48. H. Sagawa, N. van Giai and T. Suzuki, Phys. Rev. C 53, 2163 (1996) 

49. P. Navratil, B.R. Barrett and W.E. Ormand, Phys. Rev. C 56, 2542 
(1997) 

50. K. Saito and AW. Thomas, Phys. Lett. B 363, 157 (1995) 

51. B.D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986) 

52. AW. Thomas, private communication 

53. D.H. Wilkinson, Nucl. Phys. A 511, 301 (1990); Nucl. Instrum. Meth- 
ods A 335, 172, 182, 201 (1993); Z. Phys. A 348, 129 (1994) 

54. M.A.B. Beg, R.V. Budny, R. Mohapatra and A. Sirlin, Phys. Rev. Lett. 
38, 1252 (1977) 

55. J.D. Jackson, S.B. Trciman and H.W. Wyld, Phys. Rev. 106, 517 (1957) 

56. S.J. Freedman, private communication 

57. J.C. Hardy, in Proceedings of the Workshop on the Production and Use of 
Intense Radioactive Beams at the Isospin Laboratory, cd. J.D. Garrett 
(Joint Institute for Heavy Ion Research, Oak Ridge, Tennessee: Conf. 
9210121) pg. 51 



22 



